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Abstract 

We consider a generalization of the fixed and bounded trace ensembles introduced by Bronk and 
Rosenzweig up to an arbitrary polynomial potential. In the large-n limit we prove that the two are 
equivalent and that their eigenvalue distribution coincides with that of the "canonical" ensemble 
with measure exp[-nTr V(M)]. The mapping of the corresponding phase boundaries is illuminated 
in an explicit example. In the case of a Gaussian potential we are able to derive exact expressions 
for the one- and two-point correlator for finite n, having finite support. 
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1 Introduction 



Random matrix ensembles have been extensively studied since the early works of Wigner and Dyson, 
as effective mathematical reference models for the description of statistical properties in the spectra of 
complex physical systems, ranging from such diverse areas as nuclear resonances or quantum billiards 
to mesoscopic transport or quenched QCD. Even a cursory glance at some recent review monographs 
BJ , |2| , Q , |Q] , shows the impressive development of analytical tools and the variety of applications to 
physical systems reached in the past decade and a combined bibliography, although very incomplete, 
of over a thousand papers. 

Historically the matrices of the ensemble belong to one of three classes, they are real symmetric or 
complex-Hermitian or with quaternionic entries but in recent years other ensembles , like complex 
non-Hermitian or real non-symmetric matrices have been studied. To keep our paper as simple as 
possible, we restrict ourselves to complex-Hermitian matrices, although the results of this paper apply 
also to the other two traditional ensembles with minimal changes. 

A random matrix ensemble is defined by the joint probability density for the independent entries of the 
matrix. In a large number of papers, particularly those related to two-dimensional quantum gravity, 
the probability density has the form 

V(M) = L e -PTTV(M) (L1) 

where V(x) is a polynomial. Since this probability density is invariant under the similarity transfor- 
mation M = UhU* which diagonalizes the matrix M , most problems may be formulated in terms of 
the joint probability density for the eigenvalues 

P(Ai,..,A n ) = - z Al(X)e-^" v ^ (1.2) 

and may be called eigenvalue models. We shall call the probability density or the "canon- 
ical" density. In the analysis of the "large-n " limit of observables, evaluated with the "canonical" 
probability density, the method of orthogonal polynomials ^ proved to be most effective. In the 
present paper we study matrix ensembles defined by the probability density 

V S (M) = 1- 5 (A 2 -±-Tr V(M)^j (1.3) 

and the closely related probability density where the step function replaces the Dirac delta function. 
We follow the classic book by Mehta @ and call collectively these models (generalized) restricted 
trace ensembles. They are a generalization of ensembles studied long ago by Rosenzweig and Bronk 
Q where only the case V(x) = x 2 was considered. While the ensemble is still invariant under the 
unitary transformation which diagonalizes the random matrix, the method of orthogonal polynomials 
cannot be directly applied because the constraint of the delta function introduces an additional inter- 
action among the eigenvalues. Restricted trace ensembles seem to us interesting for several features 
: the interaction among eigenvalues is introduced through a constraint very similar to the non linear 
sigma model in quantum field theory, the spectral density has compact support both for finite n and 
in the "large-n " limit (unlike the usual Gaussian random model), and they relate to "canonical" 
probability densities ( |1.1|) or ( |1.2|) just in the same way as the microcanonical ensemble is related to 
the canonical ensemble in statistical mechanics. 

The effectiveness of random matrix theory is related mainly to "universal" properties of some observ- 
ables, that is independence, in the "large-n " limit, of some observables from the chosen probability 
density. The ensemble averaged density of eigenvalues p(X) = ^Tr < 5(X — M) > , in "canonical" 
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eigenvalue models (|Ll|) or (|1,2| ) is known to depend from the chosen function V(x), yet a number of 
critical exponents deduced from the spectral density were shown to be independent from the details 
of the chosen function V(x). Much earlier the density-density connected correlator p c (\,p,) 

pJX, p) = < - Trd(\ — M) — Tr5(p — M) > — < — TW(A - M) >< - Tr5(u — M) > = 

77- 77- 77 77 

= p(\,n)-p(\)p(jj) (1.4) 

was shown to have "local universality" properties , that is for |A — p\ ~ 0(l/n) and far from the 
extrema of the support of the spectral density, in the "large-n" limit. This was the basis for the use of 
random matrix theory for statistical fluctuations of observables around their mean values. Other forms 
of universality were derived more recently by several authors, including a form of "wide correlator" 
which depends on the "canonical" potential V(x) only through the extrema of the spectral density. 
The proof by Beenakker and a list of other authors is recalled in sect.l D of ref. 
In sect. 2 we exploit a scale transformation, already used by Rosenzweig in a more limited extent, to 
relate observables in restricted trace ensembles where V(x) = x 2 with the corresponding ones in the 
random Gaussian model. This allows explicit evaluations for the spectral density and the two point 
correlators for finite n. 

We then consider a generalization of the restricted trace ensembles to a generic V{x) in sect. 3. There a 
very general proof of the equivalence, in the large-n limit, of the generalized restricted trace ensembles 
with the corresponding "canonical" ones is presented. This proof is a wide generalization of the 
old result of the equivalence, in the large-n limit, of the restricted trace ensembles with the random 
Gaussian model . 

Unlike the original restricted trace ensembles, the generalized ensembles have a non trivial phase 
diagram in the large-n limit. Despite the equivalence shown in sect. 3 with "canonical" probability 
distributions, the mapping of parameters in "equivalent" models is one-to-one only in the "perturbative 
phase". We show this in detail in one example of phase diagram in sect. 4. 

Let us stress that the present paper is concerned with derivation of exact analytic results for the 
probability distributions we consider. Applications of physical interest are deferred to a future paper. 
While this paper was being written, we were informed of a poster presented by T.Nagao at StatPhys 20, 
discussing generalized fixed trace ensembles of random matrices. There the old model by Rosenzweig 
is generalized by considering a joint probability density of eigenvalues of the form 

-. n n 

V(X 1 ,..,Xn) = -Ag(A) l[\?5(A 2 -Y: A?) • (1.5) 
z i i 

This study has very little overlap with the present paper. 



2 Restricted trace ensembles at finite n. 

Let us begin by describing the most relevant features of two closely related ensembles : the fixed trace 
and the bounded trace ensembles. Let M be a n x n Hermitian matrix. The fixed trace ensemble 
corresponding to the Gaussian model is defined by the probability distribution 

V S (M) = ^-5{A 2 --TrM 2 ) , 

Z & n 

W-lTtM^fif-"^^ , (2.1) 



n 



2 



2 7T"'/ 2 

r(#) 



is the surface area of the unit 



where VM = H i=ln dMu JJ^j RedMy ImdMij , u n 2 

sphere in n 2 dimensions, and the factor ^ has been introduced in view of the large-n limit. 
Expectation values of 0(n 2 ) invariant amplitudes are trivially evaluated for every n as for instance 



< (TrM 2 ) fc > s = J VM (TrM 2 ) fc V S {M) = (nA 2 ) k 



(2.2) 



However we are interested in more general expectation values, which are functions of the distribution 
of eigenvalues. They may be evaluated from the joint probability distribution ^(Ai, .., A n ) which is 
obtained from eq.( |2.lD after integration of the unitary degrees of freedom 

1 1 n 

n(Ai,..,A„) = -A 2 (A),5(A 2 --^A 2 ), A n (A) = [] (A r - A s ) = det [Af 1 ] , 



l<r<s<n 
„2 



J i=l 



n r . 



A 2 \^2 1 n" 2 / 2 (27rf/ 



2 n 



r(f ) 



(2.3) 



Closely related to this matrix ensemble is the bounded trace ensemble. It is defined by the 
probability distribution 



6(A 2 — —Ti M 2 ) , 
Za n 



VM6(A 2 - -TrM 2 
n 



V e {M) : 

In the same way of eq.( |2.2| ), one easily finds 

< (TrM 2 ) fc > 



n (n 2 -n)/2 (Ayjn) n2 

Ur,2 o 



n 



2+k 



n 2 + 2k 



A 



2k 



(2.4) 



(2.5) 



which exhibits the usual factorization of 0(n 2 ) invariant amplitudes only in the large-n limit. In order 
to evaluate expectations which only depend on the distribution probability of the eigenvalues, one 
may use the joint probability distribution, analogous to eq. fl2.3|) : 



1 1 n 

r e (Xi,..,Xn) = -A 2 n (\)6(A 2 --J2 Af) 



f fl <iX, A»(A) <)(A* - i t A?) = (f ) n"'<°W^- f[ 
J-oo i=l n x \ z j i + ij j=1 



(2.6) 



Of course the two ensembles are related by a simple differential equation. Since 

,2 



dA 2 6K ' Za 



Vs(M)-V e {M) 



n 
2I 2 



V & (M)-V e {M) 



one easily obtains a simple relation between the two expectations for any generic observable 

2^ 2 d 



< O(M) > s = 1 + 



2 Q A 2 



< O(M) > 6 



(2.7) 
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A remarkable feature of both the "fixed trace ensemble" and the "bounded trace ensemble" is that 
the density of states p(X) has compact support for any n , finite or infinite. We here obtain the exact 
expression of the eigenvalue distribution of the fixed trace ensemble for any value of n, based on the 
known results for the Gaussian model. 

Let us first recall a few useful formulas of the Gaussian model. The partition function and the 
eigenvalue density are 

z G = Jd\ 1 ... d\ n A 2 n e -«(A?+-+^) , (2.8) 

1 r n ~ 1 
p G (X) = e~ aX2 — dX 1 ... dX n ^ At ie -«(A?+-+^-i) J] (A- Xi) 2 , (2.9) 

Z G J i=1 

where the positive parameter a is arbitrary, and for shortness we set A\ = A 2 (Ai, . . . , X n ). Both 
integrals may be computed for finite n by means of orthogonal polynomials, which in this case are the 
Her mite ones: 

Instead of evaluating the sum by means of the Christoffel-Darboux formula, it is useful for our discus- 
sion to use the expansion 

to obtain, with some simple algebra: 



k=0 



To study the integrals for the fixed trace ensemble, it is convenient to adopt the following notation. 
Let us denote by uj(n, R) the surface of the sphere in R n of radius R, and by da n the element of surface 
integration. The partition function and the eigenvalue density, for |A| < R, are: 

/oo n i n n r 

l[dX t Al(X)5(A 2 --^2^) = ^ da n A 2 n , R 2 = nA 2 , (2.13) 

-°°i=l n i=i ZU Jw(n,R) 

Ps(X) = - f f[ dX, A 2 n S(X - X n ) 6(A 2 - I f; A? ) = 
i=i i=i 

,. Tl—l 

= - \ da n ^A 2 n _ x J] (A - Xi) 2 . (2.14) 

After a change of scale, to restrict both integrals to the surface of unit radius: 

z s = ^R n2 - 2 f da n A 2 n , (2.15) 
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Pi 



r n— 1 

(A) = ^-(R 2 - X 2 )^ n2 ~ 3) / d^xA^nC-TS 

^^5 L(n-l,l) Vit 



2 - A 2 



(2.16) 



Let us first evaluate the partition function. We start from the integral expression (^^) for zg, 
and change to spherical variables with radial component r. The volume element is r n ~ 1 drda n , and 
A 2 (Ai, . . . , A n ) = r n( - n -^A 2 (X 1 /r, X n /r). Therefore we have: 



ZG 



drr n2 - l e- ar2 



I da n Al 



The surface integral is the same appearing in ( |2.15| ), and we conclude: 



Z8 = z G 



n(R^E) r 
R 2 T(4 



(2.17) 



(2.18) 



The same procedure is used in the evaluation of the eigenvalue density. In radial coordinates, the 
integral for the Gaussian density is 



ZG JO 



9 9 

drr n - 2 e- ar 



f n_1 A 

/ da n -iA n _ x J] (- - A,) 

•Mn-1,1) ~=1 r 



(2.19) 



The surface contribution is much alike the one in the expression ( [2.16 ) for ps(\). To implement this 
similarity, we introduce the expansion 

n-X 2n~2 

(2.20) 



/ da n _iA2_ 1 J[ (x - Xif = c kX k 



Since pg(A) is even in A, only the even coefficients are different from zero. The expressions for the 
densities in the two ensembles are: 



p G (X) = e 



^ 2 J_ a -i(n 2 -l)I 



n-1 



ZG 



E C2 *(Av^) 2fc r 



V - 1 



k=0 



n ,„ 2 2 , (n-l) 



W(A)= _ ( ^_V) 



n-l 

; X^/ A'-'i/.' 2 X 1 ' 

k=0 



k 



R 2 = nA 2 



(2.21) 



(2.22) 



The coefficients C2k are obtained by comparing the polynomial expression in ( 2.21| ) and the exactly 

(2£)\ 



known result (2.12): 



^(2k)\ r (a^l -it 



n 



t=k 



n 



(2.23) 



More explicitly, the spectral densities for the lowest order random matrices are 
1 



PsW 
PsW 
Ps(X) 



vr y/2A 2 - A 2 ' 
35v/3 2 A 2 
{ 3 ) 



for n = 2 



3A 4 - 2A 2 A 2 + 3A 



for n = 3 , 



32 



429vrA 14 



- -A 
4 



7 r 

2 



12^ b + 30A 2 ^ 4 - 53A 4 ^ Z + 38A 



4 ,2 



for n = 4 



2028117^/5, 2 A 2 . 



5 4 (2A) 



2:! 



375A 8 - 300A 2 ^ 6 + 4490A 4 A 4 - 5996A 6 ^ 2 + 2711A 6 



, for n = 5. 
(2.24) 
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To evaluate the spectral density for the bounded trace ensemble , for finite n , one may proceed in a 
similar way as in the Gaussian case, to obtain 

MA) = - / drr n ~ 2 / dan-iA^.! TT (- - A,) 2 (2.25) 

ze Jo Jw(n-i,i) ~r x r 



and therefore 



n— 1 / \2 \ k 



"(^^-^^S^^j , *W. (2.20) 

The same result may be obtained by inverting the differential equation ( |2.7[ ). In a similar way , it is 
possible to write the explicit expressions of the two-point correlator of restricted trace ensembles in 
terms of the known two-point correlator of the Gaussian ensemble at finite n. We obtain 

i 2n ~ 2 / 2 _ _ \ 

p G (X^) = _^ e -(A 2 +M 2 ) a -§(n 2 -4) £ Cr s{X ^y { ^y T n r s _ i , (2 .27) 

" r,s=0 \ / 

P5 (A,^) = -( J R 2 -A 2 -^)i(" 2 - 3 ) 2 <V, s AV(fl 2 -A 2 -/i 2 )-^ , (2.28) 

25 r,s=0 

where the coefficients c rjS = c S)T are defined by 

„ n-2 2n-2 

(x - y) 2 / da n _ 2 A 2 _ 2 [] " A,) 2 (y - A,) 2 = £ c r , s xV • (2-29) 

•M"" 2 ' 1 ) fe=l r, S =0 



3 Generalized restricted trace ensembles at large n 

With some generality, for an arbitrary polynomial potential V{M) = J2dkM k , where M is an Her- 
mitian n x n matrix, we define the generalized fixed trace ensemble and the generalized bounded 
trace ensemble by the two probability densities : 

P S (M) = ^5 (a 2 -- TrV(M)) , (3.1) 



V e (M) = ^-e(A 2 -- TtV(M)) , (3.2) 

where Z$ and Zq are the normalization factors, and we used the same notation of the previous section, 
where they correspond to the simplest case V(x) = x 2 . 

Both ensembles are invariant under the action of the unitary group. Therefore, when changing matrix 
the parameterization from n 2 independent matrix elements to the n real eigenvalues and the parameters 
for eigenvectors, the measures factorize into a part given by the Haar measure of SU(n) and a part 
involving only the eigenvalues. The latter provides the joint probability density of the eigenvalues, 
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the starting point for all spectral statistics. Letting <fi stand for the delta or the theta function, the 
expression for the joint probability density is: 



P^(Ai, ...,X n ) = -<p(A 2 --J2 V{Xi)) A 2 (A l5 . . . , A n ) 



(3.3) 



/n 1 n 

TJ dXi <f>(A 2 - - £ V(\i)) A 2 (A l5 . . . , A n ) . (3.4) 



i=\ i=l 



The two ensembles are obviously related by the differential equation analogous to eq.(2.7): 



V 5 {\ u ■ ■ ■ , An) = (l + Ve(Xi, ■ ■ ■ , A„) , (3.5) 

which will be used to study the properties of the bounded trace ensemble from a knowledge of the 
fixed trace one. Indeed, in this general setting, the latter is easier to evaluate in the large-n limit. 
Besides the two restricted trace ensembles, it is useful to consider also the "canonical" ensemble, with 
same potential V{M) and a parameter K: 

V(X 1 ,...,X n ) = ^e- Kn ^ v( ^A 2 (X 1 ,...,X n ) , (3.6) 

„ n 

z = Jl[dX l e- K ^ v ^A\X 1 ,...,X n ) . (3.7) 

i=l 

As it is well known, the partition function for the eigenvalues may be given the interpretation as the 
partition function of a one dimensional gas of n particles with pairwise repulsive interaction and, in the 
canonical case, subject to the external potential V(X). In the restricted trace ensembles the potential 
enters as a constraint depending on the positions of all particles. This main difference makes the 
analysis of these models difficult and interesting, especially for the issue of the universality properties 
of correlators. 

While for "canonical" models the powerful technique of orthogonal polynomials applies, giving at 
least formally and for any value of n the explicit expressions of all correlators, for the restricted trace 
ensembles we must content ourselves with the analysis of the eigenvalue density in the large n limit. 
This is easily done for the fixed trace ensemble, whose 5 constraint can be taken into account in the 
energy functional through a Lagrange multiplier. In the large-n limit, the eigenvalue configuration is 
described by a normalized density p(X), and the energy functional associated to it is 

H[p] = - J dXdnp(X)p(fj,) log |A - n\ + a (a 2 - J dXp(X)V(X)^j + (3 (l - J dXp{X)\ . (3.8) 

The saddle point configuration is the one that minimizes the energy, and is precisely the sought limit 
density p$. It solves the following equation, valid for any A inside the unknown support L of p$: 



6 Hip] 
- m 



5p(X) 



-2 J dpps(p) log |A - p\ - oV(X) - (5 . (3.9) 



A derivative in A eliminates the parameter (3 associated to the constraint of normalization, and yields 
a Cauchy-Hilbert integral equation for the limit density: 

I djJL EM = V( A) , A e L . (3.10) 
JL X — p 2 
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For any a, which is still unknown, and after having fixed a geometry for the support L (an interval, 
for example) the equation ( |3,10| ) is solved using analyticity arguments, and the extrema of L are fixed 
by the normalization condition |jl5| . Inside the family of pairs L(a) and ps{X; a) parameterized by a, 
the pair that describes the large-n limit of the fixed trace ensemble is determined by the value a = a, 
solution of the equation 

A 2 = f d\p s {X;a)V(X) . (3.11) 

JL(a) 



The number (3 of the extremal solution may be evaluated from eq. Q3.9| ) by choosing a convenient 
value of A in L. 

The density p$ so far obtained, coincides with the limit density of the canonical model ( |3.7| ), with 
parameter K = a. In the particularly simple case V(M) = M 2 , one obtains also for the restricted 
trace ensemble a limit density described by Wigner's semicircle law, with radius 2A. The energy 
functional fl3.8|) evaluated at the extremum, is 



H[p s ] = - J dXdpp s (X)p s (p)log\X- p\ , (3.12) 



where the double integral may be simplified by using the equation Q3,9| ) and the constraints: 

1 1 — 

d\dfj,ps(X)ps(n) log |A - p\ = --^aA 2 - -/3 . (3.13) 
We then obtain the large-n expression of the partition function 

zs e ~\n\^^) _ e -n*/(A*) (314) 

Since Zg = -^jZq, eq.( |3.14 ) implies 



(3.15) 



z 5 Z S n^f(A^) 

A simple check is provided by the monomial potentials V(x) = x 2k . In this simple case, the normaliza- 
tion constants zg and zs may be evaluated by a rescaling of the eigenvalues with the result |j = ^z~- 
The result ( 3.15 ) is most useful and it implies the generalization of eq. (|2.7l ) 

<0(Af)> 5 = (l + Cn ^)<0(M)> e ; Cn ^-—^—. (3.16) 



By using this equation both for the spectral density and for the density-density correlator ( |1.4| ), we 
obtain an exact equation, for any n : 

P<5, c (A,u) = p s (X,p) - p 5 {X)ps{p) = 

1 + c nQ^JPe( x ^) - ( K l + c ng^jPeW[l + Cn-g^jpe{p) = 

= pe, c (A,/i) + c n -^p 9iC {X,p) - ( c ™) 2 (^MA)^7^pe(u)^) • (3.17) 

We have not proven that the generalized restricted trace matrix ensembles have a topological expansion 
in the "large-n" limit and the factorization of invariant operators, analogous to matrix ensembles 
defined by "canonical" probability densities. The analysis of next section, where the fixed trace 
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constraint is reached as a limit of the probability density Vi(M) indicates that such properties are 
very likely. Therefore it seems reasonable to assume, as for the "canonical" probability densities, 



(3.18) 



where stands for the 5 or the 6 functions. This assumption (3.18), as well as more general assump- 
tions, together with eq.( 3.17 ) and eq.( 3.16 ) imply , in the "large-n" limit p§ ;C (A,fj,) = pg >c (X,fJ,). The 
results of this section are rather general and formal. The determination of the Lagrange multiplier a 
in eq.( 3.11| ) of course depends on the model potential V(M) in a non trivial way and on the various 
phases of the model. We provide a specific example in the next section, by the study of the potential 
V(M) = g 2 M 2 + 94 M 4 . 



4 Phase transitions 

In the previous section it was shown that, in the large-n limit, the spectral density of restricted trace 
ensembles with polynomial potential V(A, where gi are the couplings, coincides with the spectral 
density of the "canonical" ensemble with potential aV(X, gj). The scaling factor a, solution of eq. 
( p. 11 ), is actually a nonlinear function of the couplings <?j. The correspondence between the two sets 



of parameters, namely gi and agi, is one-to-one only in the perturbative phase. 
In this section we show in detail the case of the even quartic potential 

V(M) = g 2 M 2 + gi M 4 . (4.1) 

where the nonlinear relation originates different phase diagrams. To this end, we find it useful to 
consider the squared trace ensemble Vi(M) 

Vi(M) = -^exp \-l (-2nA 2 TiV{M) + (TrV(M)) 2 )] , 
Zi L V / J 

Z { = JvMexp [-Z (-2nA 2 TvV(M) + (TrF(M)) 2 )] . (4.2) 

The large-n limit of the model described by the probability distribution Vi{M) is easily found by the 
saddle point approximation. These type of models, where the exponent of the Boltzmann weight is a 
sum of different powers of traces of even powers of the random matrix was analyzed in several matrix 
models in zero and one dimension || - |14|] . The additional " trace-squared" terms were interpreted to 
provide touching interactions to the dynamical triangulated surfaces defined by the matrix potential 
TrV(M). 



For any fixed I, the model in eq.(4.2)-(4.1) is equivalent in the large-n limit to a random matrix 



ensemble with the well studied "canonical" probability distribution 
p(M) = — exp(-nTrF(M)) , 



Z = J VMexp(-nTrV(M)) , V(M) = g' 2 M 2 + ^ M 4 , (4.3) 



provided the parameters g' 2 and g'^ are suitable functions of the parameters of the model in eqs.^- 
(4.1). This may be accomplished by two equations, such as the requirement that the expectations of 
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< Tr M 2 > and < Tr M 4 > should be the same for the two probability distributions. 

On the other hand, for fixed n, in the large-/ limit, Vi(M) reproduces precisely the generalization 

of the fixed trace ensemble Vs{M) , as one sees from the following representation of the 5-function 

5{x) = lim^oo y^exp(— lx 2 ). Of course, when choosing g 4 = 0, we merely reobtain the results of the 
analysis by Bronk and Rosenzweig. 

Let us now recall the saddle point analysis for the large-n limit of the ensemble Vi{M), eqs. (|4.2j )- 
([f Since it proceeds along well known analysis, we include, for more generality the cases of the 
random matrix M belonging to the orthogonal, unitary, or symplectic ensembles, corresponding to 
the parameter [3 = 1,2 or 4. It is important to notice that, unlike the familiar quartic probability 
distribution (|4.3| ), the probability distribution ( f4.2D -([P|) is well defined for any real value of the two 
parameters g 2 , 9i- Let us begin by assuming g 2 > 0, 54 > 0, which corresponds to the perturbative 
(or one-cut) phase ; later in the section the complete phase diagram will be described. For any finite 
positive value of the parameter I, the density of eigenvalues pi(X) is the solution of the singular integral 
equation 

pj d »YZ^ = 2l(g 2 c 2 +gACA-A 2 ) V'(X) = 2g' 2 X + 4g' A A 3 (4.4) 
where the moments c/t are defined by 

c fc = J d\\ k Pl {\) (4.5) 

and g' k are the effective couplings : 

g' k = 2l(g 2 c 2 + gi c 4 - A 2 ) g k . (4.6) 

From the symmetry of the potential the support of pi(X) is expected to be one segment or two segments, 
in either case symmetric with respect to the origin. The solution of the saddle-point equation (fO|) in 
the one segment phase reads 



PiW = i~{92 + 9 4 b 2 + 2g' 4 X 2 )VW^ , (4.7) 

where the endpoint of the support [—6, b] is given by the normalization condition on the eigenvalue 
density 

f b 52 
1 = J b dXpi(X) = 21 (g 2 c 2 + 94 c 4 - A 2 ) —(2g 2 + 3b 2 g 4 ) , (4.8) 

where we have used again the g^s. 

The moments c 2 and C4 can be obtained when requiring self consistency by inserting the solution eq. 
( [4.7] ) back into the definitions ( [4.5[ ), which yields the linear system of equations 

c 2 = -21 (g 2 c 2 + g 4 c 4 - A 2 ) —{g 2 + 2b 2 g 4 ) , 

I 2/ ( 92 c 2 + 5 4C 4 - ^(4 52 + 96 2 54 ) • (4.9) 



c 4 



For a potential of higher degree we will again get a linear system of equations for the corresponding 
moments c^, fe = l, ..,m, which is due to the fact that the solution of the saddle-point equation will 
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again depend linearly on the coupling constants g' k as in eq. ( |4.7|) . Instead of solving the eqs. (|4.9| 
for C2 and c 4 we can also express them entirely in terms of the couplings with the help of eq 

b 2 (g2 + 2b 2 g 4 ) 

C2 



2(2 52 + 36 2 54 ) ' 
6 4 (4g 2 + 9b 2 g 4 ) 

C4 = 16(2 52 + 3^ 4 ) ■ (410) 

The same trick can be used to express the eigenvalue density eq. ( |4.7D only in terms of the g^, which 
reads 

" 1(A) = ^T8S50 (92 + + 294A!) ' <4U) 

where the endpoint of the support b is the root of the fourth order equation in b 2 

l(9(g 4 ) 2 b s + 20g 4 g 2 b 6 + §{(92? ~§A 2 g A )b A - 32A 2 g 2 b 2 ) =16(3 (4.12) 



which, for vanishing <? 4 and positive g 2 is asymptotic to b 2 ~ 2(A 2 + a/A 4 + j3/(2l))/g 2 . By comparing 
eq([4.4[ ) with the analogous saddle point equation for the "canonical" quartic probability distribution 
( fOj ) it is obvious that they have the same eigenvalue density, in the large-n limit, for both phases of the 
model, provided the effective coupling eq. ( |4.6| ) are precisely identified with those of the "canonical" 
distribution 

92 = 21 ( 92C2 + 54C4 " A *) 92 = b^2g 2 +3b^) 92 ' 
9', = 21 (g 2 c 2 + 54 c 4 - A 2 ) g 4 = ^J^g*) 9 * ' 

(4.13) 

Of course, the last equality on the right sides of previous equations only holds in the one cut phase. 
For simplicity, let us now proceed with (3 = 2. In terms of g' 2 and g' 4 , the equation for the support 
( [4.12 ) is the more familiar equation 3g' 4 b 4 + 2g' 2 b 2 — 4 = 0. The phase diagram of the "canonical" 



quartic model V(M), eq.(fD|), is well known: if g 2 is fixed positive , the one-cut solution (|4.7|) - ([4.8[) 
holds for any real g' 4 such that 

g',>-^(g' 2 ) 2 , (4.14) 

which is a border of existence for the model. If g 4 is fixed positive, the one-cut solution holds for any 
real g 2 such that 

92 > -2y/7i , (4-15) 
which is the line of phase transition to the symmetric two-cut solution : 

MA) = J(£>2_ A 2 )(A 2_ C 2) ( (416) 

with ends of support [-D , —C] U [C , D] being solutions of 

g' 2 +g' 4 (C 2 + D 2 ) = , g' A (D 2 - C 2 ) 2 = 4 (4.17) 
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The map between {52, 54} and {g' 2 , g'4} in this phase, may be found after the evaluation of {c 2 , C4} and 
the requirement of self-consistency , just as before. 

It is straightforward to see that the phase transition line ( 4.15| ) becomes, in the parameters of the 
model (O)-(O), the couple of lines 



94 



A 2 ± 



y/A* - 3/(2/)) , g 2 >0 , 54 <0 



(4.18) 



Therefore if A A - 3/(21) < the model (|j)-(fj) has the one-cut solution for every real value of g 2 , 
(74. In the other case A 4 — 3/(21) > the two-cut solution holds in the region of parameters bounded 



by the two curves ( 4.18} ), while the one-cut solution holds everywhere else in the plane of real values 
of #2 , 34. 

The image of the existence line ( [4.14 ) , in the space of parameters g 2 , 54 is a couple of curves : 



94 
94 



-(g 2 f(A 2 -^/A* + 7/(6l)) , g 2 >0 , <? 4 < , 



— (g 2 Y(A 2 + yjA 4 + 7/(6l)) , g 2 <0 , g 4 > . 



(4.19) 



There are two regions of the plane of real variables g 2 , 54 : the first one bounded by the first line 
( [4.19 ) and the positive axis g 2 , and the second one bounded by the second line ( 4.19j ) and the negative 
axis 52, where the equation of the support ( f4,12j ) of the one-cut solution has three possible values. The 
one-cut solution ( |4.1l| ) as function of the parameters g 2 , 54 has a first order discontinuity in these 
regions due to the cubic type instability of the solution of the eq. ( |4.12j ) with respect to the parameters 
92 1 94 ■ As usual, the lines of the first order transition are determined by comparing the evaluation 
of the free energy of the model, as functions of the different possible values of the endpoint of the 
support b. 

In the remaining part of this section, we consider the limit / — ► 00 where we obtain the distribution 
Vg(M) with the potential (4.1) explicitly. We shall denote Ck = ^ m i—>oo c k an d 6=lmn_» 00 b . Because 
of the ^-function in the distribution it will hold 



A 



92C2 + g4 C4 



(4.20) 



whereas the quantity /(52C2 + 54C4 — A ) will stay finite, as one can see from eq. ( |4.8| ). Eq. (4.20) 
is actually eq. ( |3.11| ) for the quartic potential considered in this section. Eq. (4.4) shows that the 
model with V$(M) has the same eigenvalue density of the "canonical" quartic model ( f4.3| ), provided 
ag 2 , and g 4 = ag^, where 



9' 2 



a 



lim 2 / (g 2 c 2 + g^A - A 2 



r 2 



h -(2g 2 +3b 2 g A ) 



(4.21) 



The results for the moments eqs. ( [4.10 ) and the density eq. ( [4.11 ) carry over when replacing everything 
by barred quantities. Eq. ( [4.21 ) gives the solution to eq. ( |3.11 ) and shows its dependence on the 
coupling constants of the quartic potential eq. ([O]) • 

The phase diagram for I = 00 is similar to the one previously described for finite /, with some 
simplifications. The couple of lines ( |4.18j ) which are boundaries of the two-cut phase become the line 



94 



3 (92? 
16 A 2 



92 > , 54 < 



(4.22) 
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and the negative g 4 axis. Similarly there are two regions of multiple solution for the one-cut support, 
where a first order discontinuity will occur. One is bounded by the positive part of the g 2 axis and 
the line 



9a 



7 (92? 



144 A 2 



92 > , 54 < . 



(4.23) 



The second region is the entire region 54 > and g 2 < 0. Eq.(4.20) for the endpoint b of the one-cut 
solution turns into 







9(g 4 ) 2 f + 20g 2 g 4 t + 8((g 2 ) 2 - 6g 4 A 2 )b 2 - 32g 2 A 2 



(4.24) 



The vanishing support 6 = actually provides the limiting solution /95(A) = 5(X) in the sector g 2 < 0, 
54 < 0. In other regions of parameter space the support is determined by the solution of the third 

_2 

order equation in b above. 

Let us finally extract the result for the Gaussian distribution Vs{M) with potential V(M) = g 2 M 2 



from the above formula by setting 54 = 0. Eq. ( 4.24| ) leads to 



AA 2 

92 



with the corresponding eigenvalue density from eq. (4.11) 



^ 2 ^/f^ 2 

nb 



(4.25) 



(4.26) 



this is the well known semi-circle spectral density and together with eq.( |2.7| ) it reproduces the old 
result [[?J that the spectral density of the restricted trace ensembles is equal, in the "large-n" limit, to 
the spectral density of the Gaussian ensemble. 
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